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ELECTROMAGNETIC  RADIATION  AND  HEAT  CONDUCTION 
FIELDS  NEAR  COAXIAL  CONICAL  STRUCTURES 

I.  INTRODUCTION 

The  investigation  of  the  field  induced  by  a  "vibrating"  point  or 
line  source  in  the  presence  of  conical  structures  has  been  carried  on 
for  many  years.  Most  of  the  investigations  were  concerned  with 
semi-infinite  cones  and  various  approximations  or  very  special  cases 
were  taken  to  discuss  a  finite  cone.  The  purpose  of  this  thesis  is  to 
investigate  the  finite  conical  problem  with  a  double  cone  structure  in¬ 
volved.  The  single  finite  cone  and  semi-infinite  cones  are  obtained 
as  special  cases. 

One  of  the  first  investigations  of  the  field  induced  by  a  vibrating 
point  source  in  the  presence  of  a  conical  structure  was  performed  by 
H.  S.  Carslaw  (1914).  He  divided  the  field  into  two  parts,  a  primary 
and  a  secondary,  where  one  part  satisfied  the  free- space  conditions 
and  the  second  part  satisfied  the  boundary  conditions.  The  primary 
field  obtained  involved  evaluating  c.ontour  integrals  using  residue 
theory.  The  integrands  contained  conical  functions.  The  secondary 
field  was  assumed  to  be  of  the  same  form  as  the  primary  only  con¬ 
taining  an  unknown  function  to  be  determined. 

H.  Buchholz  (1940)  examined  the  semi- infinite  ideally  conduct¬ 
ing  cone  using  the  same  procedures  as  Carslaw  used.  The  field  was 
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excited  in  this  investigation  by  an  electric  source  ring  placed  axially 
symmetric  to  the  cone's  axis. 

The  problem  did  not  have  much  appeal  until  the  advent  of  the 
space  age.  At  that  time  the  problem  of  the  finite  cone  became  of  im¬ 
portance. 

When  a  re-entry  vehicle  re-enters  the  atmosphere  at  hyper¬ 
sonic  speeds  a  double  conical  structure  of  finite  length  is  formed. 

The  interior  cone  is  the  re-entry  vehicle  itself.  The  outer  cone  is  a 
plasma  sheath  formed  from  a  surface  induced  electromagnetic  wave. 

An  electric  dipole  antennae  radiates  most  of  its  power  along 
this  sheath.  Thus  transmission  of  electromagnetic  waves  through  this 
conical  sheath  of  plasma  is  highly  attenuated.  This  attenuation  phe¬ 
nomenon  is  exhibited  by  a  marked  reduction  in  the  strength  of  the 
radio  signals  from  the  re-entry  vehicle  to  the  earth's  surface. 

There  have  been  many  investigations  of  the  field  in  the  vicinity 
of  semi-infinite  cones.  For  some  of  the  recent  investigations  one 
should  look  at  L.  B.  Felsen  (1957,  1959a,  1959b)  or  S.  Adachi  (I960). 
However,  it  is  clear  that  the  electric  dipole  results  would  be  greatly 
modified  if  a  s emi- infinite  conical  structure  were  replaced  by  a  finite 
conical  structure. 

J.  B.  Keller  (I960)  investigated  the  backs cattering  of  an  electro¬ 
magnetic  wave  incident  upon  a  perfe ctly- conducting  finite  cone.  The 
theory  was  based  upon  geometrical  theory  of  diffraction.  He  thus 
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assumes  the  results  are  "...  probably  valid  for  wavelengths  as  large 
as  cone  dimensions  or  smaller.  " 

From  1964  to  1968  several  contributions  were  made  that  are  all 
directly  related  to  this  problem. 

D.  C.  Pridemore-Brown  (1964)  investigated  the  radiation  pattern 
of  a  small-loop  antenna  situated  on  the  axis  of  a  thin  conical  sheet  of 
plasma.  He  used  the  same  basic  concept  used  by  H.  S.  Carslaw  (1914) 
and  H.  Buchholz  (1940).  He  took  as  his  secondary  field  that  induced 
by  a  surface  current  on  the  plasma  sheath.  The  eddy  current  distri¬ 
bution  had  to  be  obtained  as  the  solution  of  an  integral  equation.  By 
assuming  the  cone  was  highly  reflecting,  perturbation  techniques 
were  developed  and  applied  to  the  special  case  of  a  cone  with  apex 
angle  10°  and  the  wavelengths  of  2tt  times  the  distance  from  the 
antenna  to  the  cone  apex. 

This  was  not  too  fruitful  so  he  developed  a  different  perturba¬ 
tion  method  (D.  C.  Pridemore-Brown,  1966)  and  applied  it  to  the  inte¬ 
gral  equation  for  the  eddy  current  distribution.  He  developed  the 
equations  further  using  the  Kantorovich-Lebedev  transform  theory  to 
obtain  an  infinite  system  of  equations  involving  unknown  transforms 
instead  of  an  integral  equation.  This  process  was  first  employed  by 
A.  Leitner  and  C.  P.  Wells  (1956)  for  a  circular  disk.  Under  the  as¬ 
sumption  the  slant  height  of  the  cone  is  small  compared  to  the  wave¬ 
length,  D.  C.  Pridemore-Brown  obtained  a  solution  to  this  infinite 
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system  (1968). 

A.  Banos  et  al.  investigated  this  problem  along  the  same  lines 
as  discussed  above  (1964).  Their  solution  had  essentially  the  same 
restrictions. 

The  purpose  of  this  investigation  is  to  obtain  the  character  of 
the  field  of  electromagnetic  waves  in  the  conical  sheath  of  the  re¬ 
entering  space  vehicle  in  a  "closed  form."  The  conical  sheath  will  be 
represented  as  the  space  between  two  finite  coaxial  cones  with  coinci¬ 
dent  apexes  and  bases.  The  base  of  the  outer  cone,  and  thus  the  inner 
one  as  well,  is  taken  to  be  a  spherical  cap  of  radius  equal  to  the  slant 
height  of  the  cones.  The  base  and  outer  cone  have  a  smooth  joining 
surface.  This  geometry  approximates  the  actual  re-entry  flow  geom¬ 
etry  experienced  by  the  space  vehicle  at  hypersonic  speeds.  If  we 
assume  the  conical  shells  are  perfectly  conducting,  the  field  may  be 
investigated  by  placing  a  ring  source  of  time  harmonic  character  be¬ 
tween  the  walls  and  located  axially  symmetric  to  the  conical  axis. 

A r.  r~  '>■'/  be  seen  from  the  above,  prior  investigations  have  been 
concerned  with  determining  induced  eddy  currents,  secondary  fields, 
etc.  The  approach  taken  here  is  completely  different. 

Herein  a  very  simple  method  is  utilized.  The  Green's  function 
is  determined  by  using  the  well  known  Green's  function  representation 
as  a  series  involving  the  eigenvalues  and  eigenfunctions  of  the  prob¬ 
lem  (A.  Sommerfeld,  1964).  This  approach  leads  successfully  to  a 
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solution  for  the  finite  conical  structure  with  no  restrictions  made  upon 
the  wavelengths,  body  slant  height,  or  apex  angles. 

Chapter  II  is  concerned  with  the  derivation  of  the  eigenvalues, 
eigenfunctions,  normalization  procedure,  and  the  determination  of  the 
vector  potential  field.  The  electric  and  magnetic  components  of  the 
field  are  obtained. 

Chapter  III  examines  several  cases  of  the  solution.  Since  no 
prior  solution  has  been  obtained,  some  of  the  cases  must  resort  to 
semi- infinite  conical  forms  to  compare  with  previous  results.  In¬ 
cluded  are  the  single  finite  cone,  single  semi- infinite  cone,  semi¬ 
infinite  double  conical  wave  guide,  and  a  bi- conical  antenna. 

The  solution  of  the  problem  for  a  time  harmonic  electric  source 
ring  also  provides  the  basis  for  the  examination  of  the  heat  conduction 
inside  the  sheath.  Thus  it  was  felt  that  this  result  should  be  pre¬ 
sented  as  Chapter  IV  and  a  few  degenerate  cases  mentioned.  This 
solution  is  obtained  by  the  use  of  the  inverse  Laplace  Transform. 
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II.  ELECTRIC  SOURCE  RING  FIELD 


The  geometry  of  the  boundary  value  problem  is  determined  by 
approximating  the  actual  plasma  sheath  about  a  re-entry  vehicle. 

The  Apollo  space  vehicle  is  taken  in  this  investigation.  Thus  the  be¬ 
low  figure  depicts  the  geometry  which  must  be  used. 


The  base  of  the  two  cones  is  formed  by  a  spherical  cap  of  radius 
equal  to  the  slant  height  of  both  cones.  The  cap  is  fitted  smoothly  to 
the  outer  cone. 

The  field  between  the  cones  is  excited  by  a  time -harmonic 
electric  source  ring  of  frequency  w.  The  ring  is  placed 
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symmetrically  about  the  axis  of  the  cones. 

Let  P  =  P(r,  0,  cp)  be  any  point  of  observation  between  the  two 

cones  and  Q  =  Q(r',  0',  cp')  be  a  source  point  on  the  source  ring. 

The  source  element  at  Q(r',  0',  cp')  produces  a  vector  potential 

which  is  tangent  to  the  source  ring.  This  tangent  vector  potential  A 

has  the  three  components  A  ,  A  A  .  But  when  we  integrate  about 

r  0  cp 

the  ring  with  respect  to  cp,  and  the  ring  is  placed  symmetrically  to 

the  double  conical  structure's  axis,  the  A  component  becomes  in- 

<P 

dependent  of  cp  and  A  =  A„  =  0. 

r  0 

Maxwell's  equations  can  be  reduced  for  a  time-harmonic  elec¬ 
trical  source  to  a  form  that  says  the  vector  potential  field  u  must 
satisfy  Helmholtz's  equation 


whe  re 


A  u  +  k  u  =  0 

k  =  w/c 

c  =  speed  of  light 


Since  the  vector  potential  has  only  the  A  component,  we  may 

cp 

write 


2 

A  A  +  k  A  =  0 
cp  cp 


The  first  form  will  be  used  for  the  analysis  and  the  last  form  when 
the  field  components  are  desired. 


8 


The  electric  source  ring  induces  a  spherical  field  that  has  both 
electric  and  magnetic  components.  From  Maxwell's  equations 


E  =  v(v-A)  +  k  A,  electric  components 


H  =  V  X  A, 


magnetic  components 


Introducing  spherical  coordinates  and  noting  A=  (0,  0,  A  )  we  ob- 

<p 


E  =  E  =  0;  E  =  k  A 
r  0  <p  <p 


i  aA 

H  =  -[  — f  +  A  cote] 
r  r  90  <p  J 

i  9A 

H  =  -  -  A  -  - £ 

0  r  cp  3r 


H  =  0 
<P 


The  geometry,  physics,  and  above  equations  all  suggest  the 
boundary  conditions  are  expressed  as: 


0  =  0n,  0 .  :  E  =  E  =  0 
0  1  r  tp 


r  =  a  :  E  =E„  =  0 

<p  0 


These  suggest  that  A  =0  for  0  =  0  0  and  r  =  a.  Hence  we 

<p  0  1 

wish  to  obtain  the  first  Green's  function  for  this  field  and  geometry 
due  to  a  point  source  located  at  Q. 


The  first  Green's  function  will  be  obtained  by  using  the  well 
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known  expression  of  it  as  a  series  expansion  involving  the  eigenvalues 
and  normalized  eigenfunctions  (see,  for  instance,  A.  Sommerfeld 
1964,  p.  183): 


(3) 


G(P,  Q) 


■I 


Uj(p)Uj  (Q) 


ri  k2-k2 
(j)  J 


A  remark  on  the  notation  should  be  made.  In  general  this  sum¬ 
mation  is  usually  over  three  indices.  The  m(P)  are  the  normalized 

eigenfunctions  corresponding  to  the  eigenvalue  k.  evaluated  at  the 

❖ 

point  of  observation  P.  The  u.(Q)  is  the  complex  conjugate  of 

u.(P),  only  evaluated  at  the  source  element  Q. 

J 

Since  we  are  in  a  spherical  coordinate  system,  the  Helmholtz 
equation  can  be  examined  using  separation  of  variables.  Letting 
u(r,  0,  <p)  =  f^(r)£  (&)£  (<p)  and  substituting  into  Helmholtz's  equation 
yields  the  three  ordinary  differential  equations 


2 


(rfj)  +  (k2- 


v(v+l)  vf 
2  ’  1 


=  0 


+  cot  0 


+  (v(v+l) 


2 

_ 

2 

sin  0 


=  0 


+  p2fc 


=  0 
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where  v  and  jj.  are  separation  parameters  to  be  determined  later. 
We  have  taken  the  Q  coordinates  as  (r',^>',  0). 

The  differential  equation  for  fj(r)  is  recognizable  as  a  form 
of  Bessel's  differential  equation.  Its  solution  is  written  as 

f,  (r)  =  ~ X  C  i  (kr) 

1  2  V+V 

where  C  i  (kr)  is  any  cylindrical  function  of  order  v+j.  Since 
V+2 

we  request  u  be  finite  at  r  =  0,  we  select  the  solution  to  be 

f  (r)  =  r'"j  i  (kr ) 

1  V  +  2 

where  J  ,  (kr)  is  the  ordinary  Bessel  function  of  order  v+t  and 
argument  kr. 

The  differential  equation  for  f^(0)  raay  be  recognized  as 
Legendre's  differential  equation.  It  has  the  solution  of  any  spherical 
nanuuuu,  xuin-uon,  or  any  linear  combination  of  them.  Solutions  for 
0  real  which  are  single -valued  and  regular  are  limited  to  the  asso¬ 
ciated  Legendre  functions  of  the  first  and  second  kinds  of  degree  v 
and  order  p,  P^(cos  0)  and  Q^(cos  0)  respectively.  Thus  a  lin¬ 
ear  combination  of  these  two  functions  would  be  a  solution  for  f,,(0). 

However,  a  useful  substitution  is  made.  The  Q^(cos  0)  can 

v 

be  written  in  the  form  (Magnus  and  Oberhettinger,  1949,  p.  63) 
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_  Q^(cos  9)  =  cs  c[tt(  v+u)]P  ^  ( -  cos  0)  -  cot[rr(  v+u)]P^(cos  0) 

TT  v  v  v 

Thus 

f  (0)  =  a  P^(cos  0)  +  a  P^(-cos  0) 

2  1  v  2  v 

The  equation  for  f^(<p)  is  solved  immediately  to  be 

t  (  \ 

f3(<p)  =  e 

We  desire  the  vector  potential  field  to  be  s ingle -valued  in  cp.  Fur¬ 
thermore,  it  must  be  even  in  <p.  Thus  we  take  p  =  m  =  0,  1,  2,  .  .  . 
and 

f 3(<p)  =  cos  (rrup) 

Combining  the  respective  equations  for  f  (r),  f  (0),  f (<p),  we  ob- 

1  «  J 

tain 


(4) 


cos  (m<p) 
u  -  I 


The  substitution  process  has  involved  the  changing  of  the  (v+|)  to  v 
on  the  Bessel  function.  This,  of  course,  modifies  the  degree  of  the 
associated  Legendre  functions. 

We  are  now  ready  to  determine  the  eigenvalues  and  eigenfunc¬ 
tions.  This  is  accomplished  by  requiring  (4)  to  satisfy  the  boundary 


conditions  posed  by  (2). 
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Consider  first  the  boundary  condition  of  u  =  0  at  r  =  0. 

Since  only  J  (kr)  depends  upon  r,  u  =  0  for  all  values  of  v 
v 

such  that  J  (ka)  =  0.  If  we  denote  the  nth  positive  root  of  J  (ka)  =  0 
v  v 

by  r  ,  we  have  the  eigenvalues 
n,v 


n,  v  a 


We  need  only  take  the  positive  roots  as  J  (£)  has  the  same  roots 
for  real 

The  v-values  are  obtained  from  the  boundary  conditions  on  0. 


We  obtain  the  simultaneous  equations  of 


a  pm  ! ( c os  6  )  +  a  Pm  i(-cos  0  )  =  0 
1  v-j  0  2  v-2  u 


a,  Pm  x( co s  0  )  +  a  Pm  j(-cos  0  )  =  0 

1  v-2  1  2  v-2  f 


These  equations  have  a  non-trivial  solution  for  a^  and  a^  only  if 
the  below  determinental  equation  is  satisfied. 


Mmi(0n,0.)  = 

V-2  0  1 


Pm1(cos  0  )  Pm  i  (-cos  0  ) 

V  -W  0  V  -  p  0 


Pm1(cos  0,)  Pm1(-cos  ej 

V-V  1  V-2  1 


Since  m,  0  ,  and  0^  are  fixed  quantities,  this  is  a  transcenden¬ 


tal  equation  in  v. 
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All  the  roots  are  real  as  they  represent  eigenvalues  of  Helm¬ 
holtz's  equation.  Furthermore  since  Pm  ,!(£,)  =  Pmi(£)  we  need 

-V-2  V-2 

only  examine  the  positive  roots. 

As  the  first  positive  root  is  v  >  l/2,  r  2J  (dr)  remains  finite 

i 

for  (kr)  tending  to  zero.  This  substantiates  the  use  of  r  2J  (kr) 

as  the  solution  to  the  differential  equation  in  f^(r). 

Let  us  denote  the  fth  positive  root  of  Mmi(0  ,  0  )  =  0  by 

V  ~2  O  1 

a 

m,  £ 

Then  the  eigenvalues  for  the  boundary  value  problem  are  found 


to  be 

(6) 


T 

k  _  with  v  =  a 

n,  m,  £  a  m,  £ 


The  values  of  n,  m  and  £  are  defined  by 


m  =  0,  1,  2,  3,  .  .  .  and  n,  £  =  1,  2,  3,  .  .  . 


With  these  eigenvalues,  the  non-normalized  eigenfunctions  are  written 
as 


(7) 


_  —  ■*»  TY1 

u  =  r  2cos  (m<p)J  (  —  t  )L  (0) 

n,  m,  £  v  a  n,  v  v 


where 

(8) 


Lm(0)  =  a  Pm  i  (cos  0)  +  a0Pm1(-cos  0) 

V  1  V-2  2  V-2 


and  v 


is  determined  by  (5)  and  (6). 
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In  order  to  normalize  the  eigenfunctions  we  define  a  numerical 
quantity  N  by  the  equation 


„  6=0 
►  ZtT  r\  3-  /O  0 


\  \  \  u  r  sin  0d0drd«> 

J  n  J  n  Ja  a  n,  m,  f 
<p~  0  r  =  0  0=0 


Then  we  define  the  normalized  eigenfunctions  by 


n,  m,  f 


n,m,f  N 


We  consider  now  the  determination  of  the  normalization  factor 

2 

N.  The  triple  integral  representation  of  N  for  (7)  can  be  written 
as  the  product  of  three  integrals,  namely: 


c  ?  r  0  ?  c a  ? 

-  \  cos  ( mcp)dcp  \  [L  (0)]  sin  0d0  \  [j  (  — t  )]  rdi 
Jrs  v  J  v  a  n,  v 


P  U  2  2tt 

\  cos  (my)d<p  =  - 


where  e  is  Neumann's  number  defined  by 
m 


1,  m  =  0 


2,  m  =  1 ,  2,  3,  .  .  . 
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The  last  integral-factor  is  evaluated  using  from  Erde^yi  (1953,  Vol. 
2,  p.  90) 


I 


2  2 

uVodC  =  ^-[j  (aU-J  xl(aUJ  ,  (aC)] 

V  £  V  Vt  1  v-1 


to  obtain 


)dr 

v 


a2  ,2 


3"  [J  (t  .)- J  (t  )J  ,(t  )] 

2  v  n,v,  v+1  n,v  v-1  n,v 
2 

2  v+1  n,v,  v-1  n,  v 


From  Erdelyi  (1953,  Vol.  2,  p.  12)  this  can  be  reduced  by  the  equa¬ 
tion 

to  yield  the  final  result  of 


)dr 

v 


Tj2+i(t  ) 

2  v+1  n,  v 


The  middle  integral-factor  may  be  evaluated  also.  The  function 
m 

(0)  is  a  solution  to  Legendre's  differential  equation, 
the  formula  derived  in  Appendix  A 


Thus  using 
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I 


0 


[L^n(0)]2sin  0d0 


“[(7- v){Lm(0  n)Lm_1(0  Jsin  0  -  Lm(0  )Lm_1(0  )sin  0  } 
2v  2  vOv  0  Ovlv  1  1 


(9) 


1,en)f(x.m(e0))-i,m(e1.)f(Lm1(en)) 

2  v- 1  0  ov  v  0  v  0  dv  v-1  0 


m .  d  m 


-  Lm  (0  )~(Lm(0  ))  +  Lm(0  )/(Lm  (0  ))}] 

v-1  lov  V  1  v  lov  v-1  1 


D (m»  v,  e0>  e1 ) 


where  ~-(Lm(0  ))  is  interpreted  as  [  -^-(L.m(0  ))] 

ov  v  0  dv  v  U  v=a 


m,  £ 


For  simplicity,  the  notation  of  D(m,v,  0^,  0^)  will  be  used  for 
most  of  the  below  discussion. 

Substituting  these  integral  factors, 


N2=I^_J2  (t  )D (m,  v,  0,0) 

e  v+1  n,v  0  1 

m 


with  v  =  a  Therefore  the  normalized  eigenfunctions  are  ob- 

m,jf  6 

tained  as 


v 


n,  m,  £ 


—  y  rn 

e2  cos(m<p)J  (—  t  )L  (0) 
m  v  a  n,v  v 

- T - 

[ruD  (m,  v,  0  , 0  )]2  a  J  (r  ) 

0  1  v+1  n,  v 
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with  v  =  a  =  ith  positive  root  of  the  equation 
m,  f 


_m 


M  x(0 


0: 


V  = 


^m  / 

P  i(COS 


e  )Pm  i(-cos  0  )■ 
0  v  -  p  1 


Pm1(cos  9  )Pmi(- cos  0  ) 

V-2  1  V-2  0 


Substituting  into  (3)  and  now  making  the  coordinate  of  Q  to  be  any 
triple  (r1,  9',  <p',  )  we  get  the  first  Green's  function  for  a  finite  double 
cone  due  to  a  point  source  at  Q. 


oo  oo  oo 


_  y-  y-  XT'  E  cos[m(^-^>')]j  (rk  )L,m(0  )Ln{9' )  J  (r'h  ) 
G(P,Q)  = - ^  2 

(rr)  ira  t  =  ln  =  lm=0  2  f  2  Tn,v  2, 

Jv+l<Tn,v)D(m-V'e0-9l)[k  ‘‘"T  >  ] 


In  order  to  describe  an  axisymmetric  source  ring  we  integrate  the 
Green's  function  over  all  source  points  on  the  ring.  Integrating  and 
noting  that 


p  2tt 

\  cos  [m(?>-<p')]  cos  (<p-(p')d(p'  =  tt& 


m 


we  obtain 


•A-<p  = 


J  (r'k  )J  (rk  JL^OJL^G') 
v  n,  v  y  n,  v  v  v 


(rr')2'a2 


1  =  1  “=*  Jv+l(Tn,v,D(1-'''60-ei)[k2-<'JrL>21 


The  series  in  n  is  a  Fourier- Be s sel  series  and  in  this  case  can 
summed.  From  Erdllyi  (Vol.  2,  p.  104) 
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V  J  (r'k  )J  (rk  ) 

\  v  n,y  y  n,v 

,  J2  (T  )[k2-T2  ] 

n=l  v  +  1  n,v/l  n,v 


f  J  (kr )[ J  (ka)Y  (kr')  -  Y  (ka)J  (kr')] ;  r  <  r' 

i  V  V  V  V  V 

— n -  J 

4Jv(ka)  )  J  (kr')[j  (ka)Y  (kr)  -  Y  (ka)J  (kr)];  r  >  r' 


where  Y  (£,)  is  Neumann's  function  defined  by 
v 


Y  (£,)  =  csec  (vtt)[J  ((,)cos(  vrr)  -  J  (4)] 
v  v  -v 


It  is  more  convenient  to  introduce  the  second  Hankel  function  H'  /(t>) 
than  to  work  with  the  Neumann  function.  Using  the  definition  of 


ur  !(d  =  jv(u  -  iyv(&) 


we  obtain  for  arbitrary  parameters  f,  and  p 


U  u;  u  vb)x  (p)  -  J  (t|  )Y  (&)  =  i[j  (UH*2^)  -  J  (p)H(2)(U] 


V  V 


V  V 


Substituting  we  obtain  for  the  vector  potential  field 


oo  1  1 

V  L  <9>L  <9’> 

'■"kr'’kr) 


where 
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a  i( co s  0)  +  a  P^-4(-cos  0) 

1  v-2  2  v 

-  Jv(kr.)H<2)(lca)] 

V 

and  (kr',kr)  is  obtained  from  (kr',kr)  by  interchanging 

r  and  r' .  The  constant  D(l,v,  0^,  0  )  maybe  evaluated  from 

(9). 

The  summation  over  l  is  imbedded  in  the  determination  of 
the  fth  positive  root  ^  of  the  equation 


1/(0)  = 
v 

(12) 

,  F  (kr\  kr  )  = 
1  v 


P1  !(cOS  0JP1  i  (-COS  0 ,  )  -  p1  1  (cos  0  )PX  ,  (-cos  0  )  =  0 
V  - 2  0  V-2  1  V-2  1  V  - 2  0 


Using  the  above  results,  the  electric  and  magnetic  components  of  the 
electromagnetic  field  may  be  obtained  from  (1).  Substituting  we  ob¬ 
tain  for  the  electric  components 

E  =  kZA 

<P  <P 

Er=Ee=0 

and  for  the  magnetic  components 


•  V  e\q')  F(kr',kr) 

TT  171  \  v  1,2  V  r  /  1  v  ,  /0lT  l,n, 

H  = r — i-  /  — prr: s — -  { (v- - )  tan  (-  L  (0 

r  2(r3r>)^  ^  D(1’V’0O’V  2  2  V 


+  [v(-^--9— 2  )  +  cot(~)]a  P1  j(-cos  0)} 
sm  t)  2  2  v-j 
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oo  1  l 

v  L  (0)L  (O' ) 

H  =  — - - r  /  77F - 7T  ft~  )  ^  ?F(kr\kr)  +  2  H  (kr\  kr)} 

6  4(r3r’)2  ^  D\1>v>0o,ei)  ]’2  1)2  v 


with 


1  ,H  (kr'.kr)  =  r—  [  F  (kr'.kr)] 
1 , 2  v  3  r  1 , 2  v 


and  finally 


H  =  0 
<P 


The  field  components  are  all  real.  This  can  be  seen  from  (10) 
and  (12). 

In  what  follows,  the  vector  potential  A  will  be  used.  If  the 

<P 

magnetic  and  electric  field  components  are  desired,  the  immediately 
above  equations  can  be  used  directly. 


III.  SPECIAL  CASES 


The  vector  potential  field  representation  obtained  in  the  previous 
chapter  is  new.  Thus  a  direct  comparison  with  previous  solutions 
can  not  be  made  directly.  However,  by  taking  some  special  situa¬ 
tions,  a  few  cases  have  been  solved  in  the  past  that  allow  a  compari¬ 
son  to  be  made. 

The  majority  of  cases  that  exist  for  comparison  are  single 
semi- infinite  conical  structures.  Thus  direct  comparisons  will  be 
made  with  those  few  cases  available. 

The  cases  of  finite  single  cone  and  a  semi-infinite  double  coni¬ 
cal  structure  are  presented  not  for  comparison  purposes  but  for  com¬ 
pleteness.  Up  to  this  time  these  solutions  were  not  known. 

Before  considering  each  of  the  special  cases,  we  must  intro¬ 
duce  a  different  form  of  the  solution.  This  is  due  to  the  behavior  of 
the  Bessel  and  Neumann  functions  for  large  arguments. 

In  order  to  investigate  the  cases  of  the  semi-infinite  conical 
structures,  the  radius  "a"  must  be  taken  to  infinity.  We  must  trans¬ 
form  our  equations  to  a  monotonically  behaving  form. 

This  is  accomplished  by  transforming  the  solution  of  Helmholtz's 
equation  to  the  solution  of  the  modified  Helmholtz  equation 


22 


This  is  obviously  accomplished  by  the  transformation 


(14) 


y  =  ik  or  k  =  -iy 


Examination  of  (11)  and  (12)  show  that  substitution  of  (14)  af¬ 
fects  the  ,  F  (kr',  kr)  term  only.  From  Magnus  and  Obe rhettinge r 
1 , 2  v 

(1949,  p.  17  and  19) 


(15) 


J  (-iyC)  =  e11TVJ  (iyU  =  elv3TT//2I  (yU 
v  v  v 


tt(2)  .  .  ivrr  (1)  2i  ivir/2 

H  (-1  y£,)=  -e  H  (iyt,  =  —  e  K(yt, 

V  V  7 T  V 


Substitution  into  F  (k’r,  kr)  yields 

1  j  ^  v 


(16)  1Fv(-i\r',  - iyr )  = 


iZirv  I  (yr) 

—  rj—T  G  (\a)K  (yr1)  -  I  (yr')K  (ya)} 
tt  I  (ya)  v  v  v  v 


=  ,R  (yr1,  yr),  r  <  r' 
1  v 


vv  c  R  (yr',yr)  from  R  (yr',  yr)  by  interchanging  r 

2  v  1  v 

and  r'. 

Thus  we  have  the  vector  potential  field  with  an  imaginary  wave 
numbe  r  y 


A 

<P 


00 


=-T  y 

rr')2  L~‘ 

'  1=1 


_ ITT 

2  ( r  r' ) 


L1(0)L1(9') 
v _ y _ 

d( i ,  v,  0 0>  r 


lj2Rv(yr',yr) 


(17) 


23 


24 


(19) 


[  F  (kr*,  kr)]  =J  (kr)H(2)(kr'),  r  <  r' 

1  v  00  v  v 

[„F  (kr',  kr)]  -J  (kr')H(2)(kr),  r  >  r' 

2  v  oo  v  v 


We  are  now  ready  to  investigate  the  special  cases. 

Case  1 


The  first  case  we  will  investigate  does  not  have  a  known  exist¬ 
ing  solution  with  which  we  can  compare.  However  by  taking  a  special 
case  of  this  case  a  comparison  can  be  made.  The  two  situations  are 
treated  separately  as  they  are  both  interesting  and  useful  results  in 
their  own  right. 

We  let  0^=0.  The  finite  double -cone  structure  then  behaves 
as  a  finite  single-cone  structure.  From  (11)  we  have 


A 

<P 


oo  1  1 

V  L  <0)L  <9') 

2  i,2Fv(kr' kr) 


As  can  be  seen,  the  only  term  apparently  modified  is  the  term 
D(l,v,  0 q,  0).  However  0^  =  0  implies  only  one  boundary  is  pre¬ 
sent  and  thus  a^  =  0.  Hence,  from  (8), 


[L>»e  =o  =  PU(C0S  e) 

1  2 
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Now  a2  -  0  and  0^  =  0  provide  from  (A. 3)  in  Appendix  A  the  re¬ 


sult  of 


(20)  D(l,v,e0,  0) 


=  17  !(l-v,p!,.i(c03  e0,Pv-|(cos  %> si"  % 

+  (?+>')*P4<cos  9o>a7  (p‘-i(cos  eo))"P4(oos  0o)^(p-| (cos  1 

But  from  (13),  as  a^  =  0,  the  values  of  v  are  now  determined 


from 


(21) 


P  !  (cos  0 
v  ~2 


0 


Substituting  (21)  into  (20)  we  obtain 


D(l,v,  0 


0’ 


0)  ?i;(i+v,pi.i 


(cob  eo)-£<p!,i<cos  v 


where  —  (P  x(cos0  ))  is  understood  to  mean 

OV  V-2  0 


[£<pl 


V- 


i(cos  0  J)] 


v=a 


1,  I 


We  shall  use  this  convention  in  the  below  analysis.  But  from  Erdelyi 
(1953,  Yol.  1,  p.  161)  and  from  (21) 
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■  (^-+v)P1  j(cos  0  )  =  sin  0  P2  j_ (cos  0  )  -  (v--|)cos  9  P1  i(cos  0  ) 
Z  v  —  p  0  0  v  —  p  0  Z  U  v  -  p  0 


=  sin  0  P  x  (cos  0  ) 
0  v-2  0 


where  the  superscript  2  denotes  second  order  and  not  a  squared 
quantity.  Therefore  we  obtain  for  the  vector  potential  field 


00 


(22) 


(A  ) 


<P 


t r  i 
T 

(rr1  )2s  in  9 


Y  vU  (0,  0',  0  )  F  (kr'.kr) 
Lj  v  0  1 , 2  v 

1  =  1 


whe  re 


(23) 


V9’0'  v  = 


P1  x  ( cos  0JP1  x  (  CO  s  0’  ) 


v-2 


V-JL 


p2  x(cos  0  J  —  P1  i  ( c o s  0  ) 
v-2  0  dv  v-t  0 


r 


J  (kr) 


-{J  (ka)H(2,(kr')-J  (kr’)H(2)(ka));r<  r'  i.e.,  ,F 
J  ka  v  v  v  v  lv 


l,2Fv(kr'kr>=( 


J  ( k  r  ^ ) 

{J  (ka)H(2)(kr)-J  (kr)H(2)(ka)j  ; 
J  (ka)  v  v  v  v 


r  >  r\  i.e.,  _F 

c.  v 


k. 


The  index  H  on  the  summation  is  determined  as  the  fth  positive 
root  of 

P1  i  (cos  0  )  =  0 
V-2  0 


with  respect  to  v. 
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Now  using  a  contiguous  relation  of  the  associated  Legendre 

functions  (Erd£lyi,  1953,  Vol.  1,  p.  161,  (11)  )  the  equation  for 

U  (9  01  9  )  may  be  rewritten  as 
v  0 


uv<e’e,’0o)  = 


P  11(cos  9)P~  \  (cos  9') 
v-2 _ hzl _ 

p1  j(cos  9Jt-(P  11(cos0 

V~2  u  OV  V-J 


0 


)) 


Case  2 


This  case  is  a  further  specialization  of  the  first  case.  Here 
0^  =  0  and  a  =  oo.  We  obtain  this  case  immediately  upon  using 
(19)  in  (22).  The  vector  potential  field  becomes  for  the  semi-infinite 


single  cone 


(A  ) 

<P 


oo  Cj  (kr)H^(kr'),  r  <  r' 

T - -  Y  vUy(e,  9\  0O)<  ( 

oo  (rr')2sin0^  ^  ^  V  1  J^(kr')H^  Tkr),  r  >  r' 


t r  l 


where  U  (0,  0',  0^)  is  given  by  (23). 
v  0 

This  is  the  identical  result  obtained  by  H.  Buchholtz  (1940). 

The  electric  and  magnetic  components  of  the  vector  field  are 
obtained  by  substitution  into  (1).  The  results  compare  with  those 
found  by  F.  E.  Bornis  and  C.  H.  Papas  (1958,  p.  356-358).  Bornis 
and  Papas  discuss  this  particular  field  in  detail,  especially  its  physi¬ 


cal  characteristics. 
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Case  3 

The  two  cases  above  may  be  described  as  internal  problems, 
i.  e.  ,  the  source  ring  is  located  inside  its  prescribed  boundaries. 

The  case  of  an  external  source  ring  is  obtained  simply  by 
letting  9^  =  0  and  a  =  oo,  as  in  Case  2,  only  now  letting  0^  be 
greater  than  tt/2.  This  case  was  investigated  by  L.  B.  Felsen 
(1959a)  but  his  solution  is  in  such  a  complicated  form  it  is  not  pos¬ 
sible  to  make  a  direct  comparison.  He  has  relied  upon  the  geometri¬ 
cal  optics,  selected  contour  paths  and  use  of  asymptotics  applied  to 
the  positioning  of  the  source  ring,  observation  point,  and  conical 
structure. 


Case  4 


We  consider  now  the  double  cone  with  a  =  oo.  This  corres¬ 
ponds  to  a  semi -infinite  double  conical  wave  guide. 

The  vector  potential  field  is  now  obtained  using  (11)  and  sub¬ 
stituting  (19)  instead  of  (12).  This  process  yields 


(24)  (A  ) 


oo  2(rr')z 


“  lVilV)  r<r’ 

TT  1  \  _V _ V _  J 

J  (kr')H^(kr );  r  >  r! 
v  v 


T 


<P'~.  Lj  D(l,  v,  9n,  0,  ) 


1=1 


0’  1 


If  the  normalization  factor  D(l,v,  9^,  0  )  were  to  be  ignored, 


this  result  would  compare  favorably  with  that  which  would  be  obtained 


by  R.  F.  Harrington  if  his  solution  was  completed  (1961,  p.  281). 
Again,  the  values  of  i  are  the  fth  positive  roots  to  the  trans¬ 
cendental  equation  in  v  of 
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(13) 

M1  i(0  ,  0 

V  -  2  0 


)  =  P1  1(cos0n)P1  i(-cos  9.  )  - 

V-2  0  V-2  1 


P^  if-COS  0  )P  .(cos  0,  )  = 
v-t  0  v-p  1 


Case  5 


As  a  final  case,  the  biconical  antenna  is  examined.  This  spe¬ 
cial  case  is  obtained  by  letting  0^  =  tt  -  0^  and  a  =  co.  From  the 
geometry  it  is  apparent  we  have  again  an  external  vector  potential 
field. 

The  vector  potential  can  be  written  almost  immediately  if  one 
accepts  a  simple  substitution  of  0^  =  it  -  0^  D(l,v,  0^,  0^) 

and  the  deterministic  equation  (13)  for  determining  the  v-values. 
However,  care  must  be  taken.  Since  0  =  ir  -  0^,  we  have 


Therefore 


cos  0^  =  -cos  0^  and  sin  0^  =  sin  0^ 


M1  i  (9.,  0 .  )  =  [P1  i  (-cos  0J]2  -  [P1  i  (cos  0  )]  =  0 

V-p  0  1  V-p  1  V-p  1 


or,  the  v-values  are  the  fth  positive  roots  of  the  transcendental 


equation 
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(25) 


P1  !  (-cos  6,  )  =  P1  ,  (cos  0,  ) 
V-2  1  V-2  1 


Now  substituting  (2  5)  into  (A.  3)  of  Appendix  A  we  obtain 


=  E(v,  0] ) 


This  expression  is  derived  in  Appendix  B. 

Substituting  E(v,  0^)  into  (24)  in  place  of  D(l,v,  0  ,  0  ),  we 
obtain  the  vector  potential 


2(rr'  ¥ 


^  L1  (0  )LL  (01) 
\  __v _ v 

L  E(v,  0  :  ) 


J  (kr')H(2)(kr);  r'  <  r 
v  v 


J  (kr)H'  '(kr1);  r'  >  r 
v  v 


where  now  v  is  selected,  and  thus  £ ,  as  to  (25)  instead  of  (13). 

The  biconical  antenna  is  discussed  in  S.  Adachi  (I960).  His  re¬ 
sults  are  restricted  to  very  thin  and  very  wide-angle  cones. 

We  have  investigated  several  cases,  however  many  more  can  be 
chosen.  But  these  cases  chosen  seem  to  provide  a  spectrum  of  the 
more  interesting  situations  and/or  situations  which  provide  a  compar¬ 
ison  of  the  main  results  of  this  investigation. 
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I  i 

where  R  (r'y2,  ry2  )  is  defined  by 

1 ,  <L  v 

(27)  i 

I  i  2ie2lTiv  Iv(^r)  i  I  ii 

^(r'y2,  ry2)  =  - - - t~  {May2  )lMr'y2  )  -  Mr'y2  )fMay2  )} 

I  (y2a)  V  V  v  v 

v 


for  r  <  r'.  We  interchange  r  and  r'  to  obtain  for  r  >  r' 

l  i 

R  (y2r',  y2r).  As  is  easily  seen  this  approach  in  solving  the  heat 
conduction  problem  does  not  seem  to  be  fruitful. 


However  if  we  take  two  simple  cases  the  method  can  be  used 


successfully.  We  divide  the  class  of  problems  into  the  classes  of 


"a"  finite  and  "a"  infinite. 


"a"  Infinite 


i  JL 

In  the  case  of  a  =  oo  the  function  R  (r'y2,  ry2  )  as  ex- 

1,  2  v 

pressed  by  (27)  reduces  to 


1  1  2iei2lTV  I  I 

[  R  (r'y2,  ry2  )]  =  -  I  (ry2  )K  (r'y2  );  r  <  r' 

tv  oo  ir  v  v  ' 


A  i  2iei2TrV  A  I 

[  R  (r'y2 ,  ry2  )]  =  -  I  (r'y2)K  (ry2);  r  >  r' 

V  IT  V  V 


The  parameter  y  enters  the  vector  potential  A  only  in  these 

<P 

1  1  l 

terms.  Thus  we  need  to  obtain  ^  {R  (py2 ,  qy2  ),  y]  for  p  and  q 
arbitrary  real  positive  quantities. 


We  consider  the  transformation  of  the  form 
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4  1  11 

p  =  a2  -  p2  and  q  =  a2  +  p2 


Thus  we  now  require  p  >  q.  Solving  for  a  and  (3  we  obtain 


a  =  |(p+q)2  >  0 


and  (3  =  ^(p-q)2  >  0, 


as  p  and  q  are  both  real  positive  quantities.  Since  a  and  (3 
are  both  positive,  we  obtain  from  Erdelyi  (1954,  p.  284,  (56)  ) 


£>  ^  (py2  )Kv(qv2 );  y) 

t 

=  oC  [(cn/)2 -(py)2  ]K  [(ay)2  +  (Py)2  1 !  y) 


t- "  Kt 


I  <f^> 

V  2Kt 
2Kt 


(q+P) 

2Kt 


t~*  Kt 

.  2  2 

(p  +q  ) 


_L  i  (_ESL  \  e  4Kt 

2Kt  V  4k t  ’ 


We  note  the  resulting  inverse  Laplace  transform  is  symmetric 
in  p  and  q.  Thus  the  first  Green's  function  in  a  semi-infinite 
double  -  conical  structure  may  be  written  as 


(28) 

G(P,  Q,  t)  = 


oo  1  1 

L  (0 )L  (0') 

-TT  1  \  V  V 


2(rr')2  D(1,V)  VV 


2ieZirlV  rr' 

£iE -  i  (_)  e 

2Kt  V4Kt' 


,  2  A 

(r  +r'  ) 
4Kt 


tr  e 


/  2  A 

(r  +r'  ) 
4k  t 


2Kt(rr')2 


oo  1  1 

L  (0 )L  (0)  _  . 

v  v  2iriv  .  rr' 

D  ( 1 ,  v,  0  ,  0  )  v(4^t) 


=  1 


O’  1 
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where  the  index  £  is  connected  with  the  v-values  as  the  £th 
positive  root  of  the  transcendental  equation  in  v  of 

(13)  p1  i  (cos  ejp1  ji-cos  ej  .  p1  !  ( -  cos  ojp1  j  (cos  ej  =  o 

V-2  U  V-2  1  V  -  ^  U  V-2  1 

The  quantities  I_/(0)  and  D(l,v,  0  ,  0^  )  are  determined,  as  be¬ 
fore,  from  (12)  and  (9).  We  see  that  the  solution  has  the  typical  ex¬ 
ponential  decay  behavior. 

The  cases  examined  in  Chapter  III  for  the  semi- infinite  conical 
structures  may  be  re-examined  again  only  for  the  heat  conduction 
instead  of  the  vector  potential  due  to  a  time  harmonic  electric  source 
ring.  The  only  difference  would  lie  in  use  of  (28)  instead  of  (1  1  ). 

The  specific  case  of  a  single  semi- infinite  cone  (Case  2  of 
Chapter  III)  was  investigated  by  R.  Muki  and  E.  Sternberg  (I960). 
They  approached  solving  the  problem  by  using  the  Mellin  transform 
and  integration  by  parts  to  reduce  the  partial  differential  equation  to 
an  ordinary  differential  equation  of  the  Sturm  -  Louiville  type  in  terms 
of  the  transformation.  The  solution  obtained  is  in  the  form  of  evalu¬ 
ating  the  inverse  transformation.  This  leads  to  an  improper  integral 
whose  integrand  contains  the  ratio  of  two  conical  functions  whose  de¬ 
gree  involves  the  variable  of  integration.  The  solution  must  be  ob¬ 
tained  by  using  numerical  integration  techniques  of  a  result  obtained 
from  numerical  integrations  of  the  conical  functions.  This  is  quite 
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cum  be  rsome. 

Application  of  the  same  processes  used  in  Chapter  III  provide 
us  with  the  solution  in  closed  form. 

Letting  0^=0,  the  below  equation  is  valid: 


L1(0)L1(0') 
v _ v _ 

Dd.v.e/j) 


P“\ (cos  0  )P~  \  (cos  0') 

V-2 _ v-j _ 

P  !  (cos  0  [P~  \  (cos 

V  -  2  0  dv  V-2 


2  v 


sin  0 


0 


u  (0,  0',  0  ) 

s  m  0  ^  v  0 


Substituting  into  (28)  we  have  the  representation  of  the  heat  conduc¬ 
tion  in  a  semi-infinite  single  cone  expressed  by 


G(P,  Q,  t) 


,  2  2 
(r  +r'  ) 

4Kt 
tt  e 

- j - 


Kt(rr'  )2sin  0 


l*V*’ 

H  =  1 


where  under  these  conditions  the  correlation  between  f  and  v 
is  given  by  the  i  being  the  £th  positive  root  of  the  transcendental 
equation  in  v  of 


(cos  0 


0 


)  =  0 


The  other  cases  examined  in  Chapter  III  for  a  -  oo  can  be  ex¬ 


pressed  similarly  for  the  heat  conduction  problem. 
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"a"  F inite 

The  problem  of  obtaining  an  inverse  Laplace  transform  of  (27) 
was  mentioned  earlier.  But  upon  investigation  of  where  (27)  originates 
shows  that  we  obtain  it  by  transforming  the  summed  series  in  n  to 
a  function  of  y  instead  of  k.  Returning  to  this  point  we  have 


oo  oo  1  1 

L  (0 )L  (0’)  J  ( r'k  )J  (rk  ) 
v _ v _  v  n,  v  v  n,  v 


%  L  D(  1,  v,  0  ,  6  )  2  ,r, 2  ,2  , 

•rra  rr1  p  „  ,  ,  0  1  J  ,(t  )lk  -k  I 

£=\  n=l  v+1  n, v /L  n,vJ 


whe  re 


k  =  t  /a 
n, v  n,v 


and  v  and  £  are  related  by  £  being  the  fth  positive  root 
of  (13). 

Substitution  of  k  =  -iy  and  using  (28)  gives  the  first  Green's 
function  for  the  heat  conduction  between  two  coaxial  finite  conical 
structures  as 


G(P,Q,t)  = 


+2 


00  00  1  1 

L  (0 )L  (0')  J  (r'k  )J  (rk  ) 
v  v  v  n,v  v  n,v 


Tra2(rr') 


L  I 


f=l  n=l 


J  ,,(t  ) 

v  +  1  n,v 


x  £ 


-1 


Y 


Kt 


From  Erdelyi  (1954,  P.  229,  (1)) 
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t 


-  1 


t“ *■  Kt 


e 


-  Ktk 


2 

n,  v 


Thus 


G(P,Q,t)  = 


•rra  (rr1)2 


,  .  - Ktk 

oo  oo  1  1  n  v 

L  (0)L  (0*)J  ( r'k  )J  (rk  )e 
v  v  v  n,  v  v  n,  v 

,  ,  Dd.v.e  e  )j2  (T  ) 

=  1  n=l  0  1  v+1  n,  v 


where  the  index  H  is  the  £th  positive  root  of  the  transcendental 
equation  in  v  expressed  by  (13). 

The  special  cases  for  finite  "a"  examined  in  Chapter  III  can  be 
treated  easily  using  the  above  equation. 
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APPENDIX  A 


The  normalization  factor  for  determining  the  normalized  eigen¬ 
functions  involves  the  evaluation  of  a  definite  integral  with  the  pro¬ 
ducts  of  associated  Legendre  functions  of  the  same  order  and  degree. 
These  relationships  are  not  available  in  the  literature  and  thus  are 
derived  in  this  appendix. 

The  Legendre  differential  equation  is  written  as 


(A.  1) 


(1-x2)  — j  -  2x  ~  [v(v+l)- 
.  £  dx 

dx 


and  has  the  solutions  of  PP(x),  Q  or  any  linear  combination  of 

v  v 

these  two  linearly  independent  transcendental  functions  over  the  real 


interval  (-1,  1).  We  take  p  and  v  as  completely  arbitrary 
paramete  rs . 

Let  wP(x)  and  vP(x)  be  any  solutions  to  the  differential 
Y  o' 

equation  (A.  1),  with  v  and  p  replaced  by  cr,  p  respectively 

for  vP(x).  Then 
cr 


dL11  dw^ 

( 1  -x^)  - t—  -  2x  V  +  [v(v+  l)-n2(l-x2)  =  0 

-  c  dx  v 

dx 


d2vP  dvP 

(1-X2)  - y  -  Zx  -  +  [(r(or  +  l)-p2(l-X2)  ^  ]  VP  -  0 

_  c*  ax  cr 

dx 


42 


Multiplying  the  first  equation  by  and  the  second  by  w ^  and 

then  subtracting  the  resulting  equations  gives  the  single  equation 


dw 


dw 


cr  dx  dx  v  dx  dx 


+  [v(v+ 1 )- c(c  +  l )  +  ( 1 -x^)  ^  (p^-p^)]w^v^  =0 

v  c 


This  can  be  rewritten  as 


dw^  dv^ 

d  /,  2  r  p  v  p  ^ 

~  (1-x  (vr  — ; -  - — J 

dx  (r  dx  v  dx 


[  (v  -  c)(v+  c  +  l)+(p^-p^)(l-x^)  *] 

V  O' 


Integrating  both  sides  we  obtain,  for  -1<  a<b<  1, 


^  [(v-c)(v+  c+l)  +  (p2-p2)(l-x2)  *  ]w^  (x)v^  (x)dx 

a 

dv*3  dw^ 

/,  £\r  l1  o'  P  v  t 

=  (1-x  w 1  -r—  -  — —  \ 

v  dx  a  dx 


[x(v-or)w^vf3  +  (cr+pjw^v^  -(v+p)vpw^  ] 
VC  V  or-  1  CV-1 


p„p 


rP 


using  the  recurrence  formula  from  Erdelyi  (1953,  Vol.  1,  p.  161). 

Now  letting  p  =  p  and  v*3  =  w^3  we  obtain 

0"  0" 
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b 

a 

If  v  =  or,  the  second  term  has  an  indeterminant  form.  Thus,  for 
o-  /  -  1/2,  we  have  from  1 ’Hospital's  rule 


(A.  2) 


P  P  A 
w  w  dx 


xw^w^  (cr+pjw^w*3  -  (v+t 
v  cr  v  cr- 1 

(v+cr  +  1)  (v- <r)(v+cr+l ) 


(v  +  p)w^  w^ 
v- 1  cr 


lim 
v-"  cr 


(cr+p)w^w^  -(v+p)w^  , 

v  cr- 1  v- 1  cr 

(v  -  cr)(v+cr+ 1 ) 


.  .  3w^  9w^  w^w^ 

l?+£l.  rwP  _ Z  _  WP _ Zzl)  .  17  .  °r- 1 

(2cr+l)  cr-1  3cr  cr  3cr  2cr+l 


Combining  this  limit  with  (A.  2)  for  v  =  cr  we  obtain 


where  w^(x)  is  any  solution  of  Legendre's  differential  equation  on 
cr 

the  cut  -  1  <  x  <  1 . 

Chapter  II  required  the  evaluation  of  an  integral  of  this  form. 
Since  (8)  is  a  solution  of  Legendre's  differential  equation,  we  can  sub¬ 
stitute  into  the  above  result  to  obtain 
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n  0  1 

[L W 


sin  9  d9 


i  f  i  ?  ll  1  r  1  ^(9)  ,  9L1  (9) 

~  /  [L1  (9  )]  cos  9-L  (Q)L  (9)+(-  +  v)L  (9)— —  -L(9)-~— 
2v  \  v  v  v-1  2  v -1  ov  v  ov 


Now 


cos  9[L1(9)]2  -  L1  (9  JL1  ,(9)  =  (i-v)  sin  9L1(9)L°(9) 
v  vv-1  2  vv 


Thus 


I 


sin  9  d9 


~  <(  j-v)[L](9  )L°(9  )  sin  9  -  LJ(9  )L°(9  )  sin  9  ] 

2v  i  2  v9v9  9  v  1  v  1  1 


i+v)  [^“V 


,  aiAe I  )  aP  (0 )  ad(e)  ad  fe,) 

v-W—  -  W5F-  -Lv.i(e>>—  +  Lv<e,>-f7“ 


From  (8)  we  obtain  the  desired  integral  solution.  By  making  the  sub¬ 
stitutions  of 


P1  !  (cos  9)  =  P1  (9) 
V-2  V- 


we  obtain  the  formula  expressed  as  (9)  only  in  its  expanded  form 
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(A.  3) 


-2vD(i, v,  0Q,  el ) 


=  -  2v 


p  ®  1  7  P  ®  1  12 

\  [L  (0)]  sin  0d0  =  -2v  \  [aP(0)+aP  (-0)]sin0d0 

Jn  v  Ja  1  M-  2  p 


{ 1  +p) 


=  (0o>sin  0O-P|iL(ei  >P^(61 )sin  91] 

r ,  apV)  ,  ap1  ,(en)  ,  apta)  .  ap1  ,(e. n"\ 
+  <1+{P,1(9^  -p,4(V-^+ 

+  a  a  J(-n)[pV)P  (-0J+P1  (-en)P  (0J-P1  (0  )P  (-0  )-P2  (-0  )P  (0, )] 

1  2\  p  0  p  0  p  0  p  0  pip  1  p  1  p  1 

r  i  apVej  ,  dp\e)  apV© )  .  ap1  (0 ) 

,  3P1  ,(-01  ,  3P1  (01  SP1  (-6  )  3P1  (0  )-|'l 

-  P^-^—  — +V-ei)^ — J/ 

+  »lii‘<Pi<-9o)pl.«-eo>  si“  eo  -  p(l<-ei>ptl<-ei)  si"  V 

r- 1  apl(-0n)  1  ap1  (-0  )  dp\.e  )  ap^-ej 

+  (1+ll)[Pp-l,"eo)_3v  "Pn("eo)-^  -p,-1<-ei>^-+y-e.>-V- 

This  expanded  form  is  quite  useful  for  examining  degenerate  cases. 


fi  d 

■>|_vw- 
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APPENDIX  B 

The  biconical  antenna  case  for  the  comparison  of  the  main  re¬ 
sults  of  this  investigation  with  known  results  requires  the  evaluation 
of  D(l,v,  ir-0j,  0^),  which  was  redefined  as  E(v,  9^). 

Substituting  0^  =  tt  -  0^  provides  the  equations  of 

(B.  1)  cos  9  =  -  cos  9^ 

sin  0^  =  sin  0 

Thus  from  these  equations,  (13)  is  rewritten  as 

[p1  i  (-cos  0  )]2  -  [P1  x(cos  0  )]2  =  0 
V - g  1  V-2  1 

Thus  the  roots  are  obtained  from 

(B.  2)  P1  i  (-cos  0  )  =  P1  i  ( co s  9  ) 

V-2  1  V-2  1 

Using  the  notation  introduced  in  the  previous  appendix,  i.  e.  , 

P1  ,(003  0)  =  P1  (0), 
v-2  H- 

substitution  of  (B.  1)  into  (A.  3)  provides  the  equation 
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-  2vE(  v,  0  )  =  afA-nJtP^-e  )P  (-e.)sin  0  -PJ(0  )P  (0  )sin  0  ] 

1  1  \  (JLiLil  1  (1  1  11  1  1 


+  ( 1  +  fx) 


[p1  ,(-e 

[_  [L-l  1 


dP1  (-0.  )  ,  3P1  ,(-0  ) 

) - f -  -  P  (-0  ) - ^ - - 

a  v  jj.  1  o  v 


3p‘(0  )  3P1  (8  )0 

-v»(v-£r- +  vv-  ■>;  / 

r  J 

+  a,/  <-H )[  p1  (0  )P  (0  )sin  0  -p\-8  )P  (-0  )sin  0  ] 
2  t  (i  1  |jl  1  1  (i  1  p  1  1 


+  (1  +  p) 


Pl  T(01 
|JL  -1  1 


a  i  ,  3P  ,  (©. ) 

)^-pl(0  )-pl(0  )  %■■■ 

OV  [i  1  }ji  1  0  V 


9P  (-0  ) 

Pl  ,  (-0!  )^-pl(-01  )  -pl(-0!  >  - - L 

(i-l  lov|i  1  |i  1  3v 


Substituting  (B.  2)  and  combining  terms  one  obtains 


,  2  2^ 
az-a  1 


E(v,6  )  =  — - - /  (-n)P  (0  (sine  (P  (.8  )-P  (0,)) 

1  iv  \  ii  1  1  u  1  n  1 


p1  (0  l-P1  (-0  ) 

(i-r  (i -r  r 


1  ->  a  PDn 

+  (1+nMP1  (0  ))  —  -tp-i - ) 

v  1  8vL  p1  (e  ) 

(1  1 


Now  from  (B.  2)  and  Erdelyi  (1953,  Vol.  1,  p.  148) 


3P  (cos  0  ) 


BP  (-cos  0  ) 


1  (cos  0  )  =  -sin  0,  — rr - =  sin  0,  — r - :  —  =  P1 


Pp(cos  V  =  "sin  91  -Tc^”0“-) 


'  -PT - HT"  =  p  (-cos  ei  ) 

1  B (cos  0  )  |i  1 


Thus 
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a 

3(cos  ) 


[P  (cos  0,  )  +  P  (-cos  0 
V-  1  V- 


This  implies  either 


P  (cos  0,  )  +  P  (-cos  9_ )  =  constant 

(JL  1  |1  1 


for  all  p  for  any  given  or 


P  (cos  0,)  +  P  (-cos  0,)  =  0 

pip  1 


Since  the  first  case  is  absurd,  thus  we  obtain 


P  (cos  0  )  =  -P  (-cos  0.) 

pip  1 


Substituting  we  obtain  the  desired  result  of 


E(v,01)  = 


,  2  2* 
2  v 


2(-p)P2  (9,  )P  (0  )sin  0, 

p  1  p  1  1 


P 


1 

p-1 


(v-p,! 


< 


pi  = 0 


(1  +  viKP1  (6  ))2 
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cone,  and  biconical  antenna.  The  heat  conduction  problem  for  the  identical 
geometry  is  also  solved.  The  solution  is  obtained  by  employing  the  theory  of 
Laplace  transformations  on  the  Green's  functions  for  the  electric  source  ring. 
Two  cases,  finite  and  semi-infinite  slant  heights,  present  two  forms  of  solu¬ 
tions.  Finite  slant  height  yields  a  double  infinite  series  and  semi- infinite  slant 
height  a  single  infinite  series. 
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